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We study the linear response and the dynamic structure factor of weakly interacting Bose gases at 
low temperatures. Going beyond lowest order in the weak coupling expansion allows us to determine 
the contribution of the thermal and quantum depletion of the condensate to the dynamic structure 
factor. We find that the quantum depletion produces a pronounced peak in the dynamic structure 
factor, which allows for its detection via a spectroscopic analysis. 



The main characteristic of bosonic quantum gases is 
the formation of a condensate at low temperatures as 
observed in 4 He and atomic gases 0, In superfiuid 
4 He, strong interactions produce a large depletion of the 
zero-momentum state and the condensate involves only 
a small fraction of the particle density; the condensate 
was observed experimentally by measuring the dynamic 
structure factor with neutron scattering j3, 0] ■ In turn, 
interactions are weak in atomic gases and the conden- 
sate involves nearly the full particle density, while the 
quantum depletion is small. However, recent progress 
in tuning atomic gases allows to increase the interaction 
strength via Feshbach resonances Q or optical lattices 
|3- It is this strongly correlated regime with large 
quantum depletion which is currently attracting a lot of 
experimental and theoretical interest. In this letter, we 
study the contribution of the quantum depletion to the 
dynamic structure factor, and propose that, in analogy 
to superfiuid 4 He, spectroscopic studies provide a power- 
ful tool for characterizing the structure of atomic gases 
and allow for the detection of the quantum depletion. 

The dynamic structure factor of an ideal bosonic gas 
exhibits a delta function peak S(u>,\i) = n 5{oj — ek), 
thus probing excitations from the condensate with the 
energy ek = ft 2 k 2 /2m, see Fig. ^a). Previous studies 
of the dynamic structure factor in atomic Bose gases ac- 
counted for the modification of such single-particle exci- 
tations via interaction and finite trapping: while the in- 
teraction modifies the quasi-particle excitation spectrum 
and the weight of the delta function peak 0, EJ , the 
finite size of the trap induces a broadening of this peak 
[ToL lllj. In this letter, we determine the linear response 
function and the dynamic structure factor of weakly in- 
teracting Bose gases within a microscopic analysis going 



to first order in the interaction parameter 



(here, 



n is the particle density and a denotes the scattering 
length). We find, that the structure factor contains two 
contributions: The first accounts for excitations of quasi- 
particles from the condensate, see Fig.^a), and involves 
the condensate density no. The second term is due to 
excitations of two quasi-particles and accounts for the 
quantum depletion n D — 1 — no of the condensate, see 
Fig. ^b); particle- hole excitations [see Fig. QIc)] only 
contribute to the structure factor at finite temperatures. 



The quantum depletion n D is well understood within the 
weak coupling Bogoliubov theory and accounts for the 
interaction-induced expulsion of particles from the zero 
momentum state. Although the extension of this theory 
to finite traps is successful in explaining many phenom- 
ena observed in cold atomic gases, the detection of the 
quantum depletion n D itself has not been achieved so far. 
We show below that the quantum depletion produces a 
pronounced peak in the dynamic structure factor allow- 
ing for its observation via a spectroscopic analysis. 




FIG. 1: Diagrams representing the excitations created by 
the external drive SV (wiggly line), (a) Creation of a quasi- 
particle (straight line) from the condensate no (zigzag line), 
(b) Creation of two quasi-particles probing the quantum de- 
pletion, (c) Creation of a particle-hole pair probing the ther- 
mal depletion. 

The microscopic Hamiltonian of interacting bosons 
takes the form 



H 



dx t/j + (x) 



h 2 A 

2m 



V>(x) 



(1) 



dxdy U(x- y)'0 + (x)V' + (y)V'(y)V'(x) 



with ip + (x) and ip(y) the bosonic field operators and 
the interaction potential U(x). In the low-density limit, 
the microscopic interaction potential U(x) is replaced 
by the s-wave scattering potential U(x) = gS(x) with 
g = Airh 2 a/m [II El 03. The expansion parameter in 
this weakly interacting Bose gas is V no 3 . 

The response Sp(t, x) of the bosonic density relates to a 
small external potential 5V(t, x) via the integral relation 

Sf(t, x) = y" dt'dx' x (t - (', x - x')SV((', x') (2) 

with the linear response function |15| (p = ^> + i/>) 

X(f,x)=-ie(f)<[p(t,x),p(0,0)]>. (3) 



2 



Here, (. . .) denotes the quantum statistical average at 
fixed temperature T and chemical potential fi. The re- 
sponse function ^(t, x) is conveniently calculated via its 
relation to the density-density correlation function in 
imaginary time D(r,x) 0, 



£>(r,x) 



<Tp(r,x)p(0,0))-(p)<p) (4) 



with T the time ordering operator. Then, the response 
x(w, x) is the proper analytic continuation of the density- 
density correlation function D(p, s , x), which respects the 
retarded character of x(i,x). Here, 57 s = 2irTs denote 
the Matsubara frequencies with s 6 Z. 



°o = 





a! 



D' = 



FIG. 2: Diagrams contributing to the response function D 
to first order in V na 3 . Do collects diagrams involving the 
condensate density no, while the diagrams in D 1 probe the 
depletion of the condensate. The thick solid line denotes the 
normal and anomalous Green's functions Q, T + , and T within 
Bogoliubov approximation. Zigzag lines account for the cre- 
ation £ + and annihilation £ of particles from the condensate. 

In the following, we apply standard quantum field the- 
ory methods for interacting bosons at low temperatures 
Ha, llM UJj • We separate the bosonic operators £ + and £ 
accounting for creation and annihilation of the zero mo- 
mentum state from the field operator, 

V + (t, x) = e + ^' + (r, x), 1>( T , x) = £ + V'(r, x). (5) 

Within each diagrammatic expression the operators £ 
and £ + are considered as c-numbers contributing a factor 
y/no; no denotes the condensate fraction and has to be 
determined self-consistently. The difference n D = n — no 
between the averaged particle density n = (p) and the 
condensate fraction uq is due to the thermal and quantum 
depletion of the condensate. Here, we focus on the low 
density limit and expand the density-density correlation 
function J3J in the small parameter V na 3 ; the relevant 
diagrams are shown in Fig. [21 All additional diagrams 
involve a vertex operator T and contribute to higher or- 
der terms in the expansion parameter y/ na 3 ; they will be 
dropped in the following. 

The diagrammatic representation of the leading con- 
tribution Do to the density-density correlation function 
D is shown in Fig. and takes the form 

D a (n s ,k) = -n (g'(fi s ,k) + g'{-n s ,-k) (6) 

+jr+(fi S) k)+jr(n S) k)). 

The normal Green's function Q' and the anomalous 
Green's functions T + and T are defined as 



0'(r,x) 



T^'(t, x)V>' + (0,0) 



^+(r,x) = -i(T^f(r,x)/+(0,0)), (7) 
no \ I 

^(r,x) = - — (t£+£ + ^(t,x)^(0,0)\. 

nn \ / 



Note, that Q' is related to the depletion of the conden- 
sate, n D — Q'(t — ► + ). Within first order in the ex- 
pansion parameter \J na 3 , the Green's functions take the 
Bogoliubov form 0, 



if] s + e k + jU 
(zft s ) 2 -£ 2 ' 



(i n s ) 2 - ^ 



(8) 
(9) 



with the excitation spectrum E k = + 2e k Ji and the 
chemical potential fi = gn. The analytic continuation of 
Do (p, s , k) provides the contribution 



Xo(w,k) 



2e k n 



(Hlu + isy 



El 



(10) 



to the response. The structure factor S(u>, k) is related 
to the imaginary part of the response function, Im% = 
-tt[1 - exp(-0hJ)]S, and using (JTUJ) we find [fi = l/T) 



5o(w,k) = 



n e k 
e^-lE^ 



5(hLo-E k ) + 5(hu + E k ) , (11) 



accounting for the excitation of single quasi-particles 
from the condensate, the process illustrated in Fig.^a). 
Furthermore, So(w, k) nearly exhausts the /-sum rule 



/DO 
w So(ijJ, q) 
-oo 



n ft 2 q 2 
2m 



(12) 



the missing part n D fi, 2 q 2 /2m derives from the structure 
of the thermal and quantum depletion of the condensate. 
With n D oc ny/na 3 , its effect is usually ignored to lowest 



order in the expansion 



but with increas- 



ing interaction its contribution plays a substantial role. 

The first diagram in the correction D 1 — D' + D' 
see Fig. [21 takes the form 



(2)) 



i»; 1) (n„k)=-T^; d q 



iVtt + e^ + u iVt s+t + e k+q + (j, 



s+t) 



E l+q. 



The term cx /i 2 diverges logarithmically for k — * 0; this 
divergence is removed in the full perturbation theory 
[l2l flq | and results in a higher order contribution of or- 




dcr 



; ln(r 



Similarly, the second diagram D' in 



Fig- contributes a correction of order na 3 ln(na 3 ). In 
the following analysis, we drop these terms and account 



for corrections of order 



Summation over Masub- 



ara frequencies fl t and analytic continuation provides us 
with 



3 



dq 1 I Eg + e q + /z 
(27r) 3 2i; q | ef* E *-l 



£q+k + M 



[r^ + 2 ^ + £; q ] 2 -£;2 +k 



-/kj + Eq + e q _ k + /j, 
[huj + id - Eq} 2 - El_ k 



(13) 



-E a 



tko - Eg + e q+k + /I -huj - Eg + e q _ k + fx 



1 - e-P E « [[hiu + iS- Eq} 2 - £2 +k [fa + iS + Eq} 2 - ££_ k _ 



The poles account for the creation of two particles out 
of the condensate with energy hut = E^ + -Ek+q [see 
Fig. db)] and the scattering of thermally excited quasi- 
particles fku — E^ — -E k +q [see Pig. Dfc)]. This creation 
of particle hole-pairs dominates in the high temperature 
limit T ^S> [i- Then, the response derives from Eq. Ijl3|l 
taking the limit /j, — ► and reduces to the bosonic Lind- 
hard function 



X>,k) 



dq f [ej - / [e q+k ] 



(14) 



(2n) 3 HlO + 6q - e q+k + iS 

with / = l/[exp(/3e)— 1] the bosonic distribution function 
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FIG. 3: Density plot of the dimensionless structure factor 
hu)S' /2/i (|15pi for an interacting gas at zero temperature; the 
momentum fe M is measured in units y/rnjl/h. The white line 
denotes the energy threshold for the two-particle creation. At 
low momenta |k| < 2^/mp/h, the peak in the structure factor 
appears at Two ~ 2/j,, while at high momenta |k| > 2^/mfi/h 
this peak follows the Bogoliubov energy huu = E^. 

Here, we are interested in the low temperature limit 
T -C fj, of the response, where thermally excited quasi- 
particles are quenched and the response function l|13(l 
accounts for the quantum depletion alone. We set T = 0, 
and obtain the correction to the dynamic structure factor 



S"(w,k) = n D 



h 2 k 2 1 



2m 4/i 2 



5" 



4/? 2/j, 



(15) 



with n D = (8/3 y^) n V na 3 the quantum depletion 
within Bogoliubov theory [l3T |. The dimensionless func- 
tion S' (a, x) vanishes below the two-particle excitation 



threshold x 2 < a 2 + 2a, while for x 2 > a 2 + 2a we obtain 

1/2 



S' (a, x) = 



8V2 a 3 / 2 



de 



1 - 



1 



(a + l+e) -x 2 



where £~ = — 1 — a + \/l + x 2 . The upper integration 
limit e + is determined by E\ < e + = £2 < £3 with £1,2,3 
the solutions of the cubic equation 



[a + e + 1) - 1 + 4a£ 



x 2 + 4ae [a + e+iy = 0. 



Combining the results l|ll(l and (|15|l provides the dy- 
namic structure factor S(u>,]i.) — 5o(w,k) + S"(o;,k) 
within first order in the expansion parameter %/ na 3 . Us- 



ing the relation dx xS'(a, 



1, the /-sum rule is 



satisfied as required in a consistent calculation to order 



Vna 3 . The dimensionless function Hu>S'/2/j, is shown in 
Fig- El At low momenta \k\ < ^/m/i/h, the main contri- 
bution appears at high energies hw ss 2fj, and accounts 
for the excitation of two quasi-particles with opposite 
momenta, see Fig.^b). This peak at low momenta is a 
characteristic of the quantum depletion in an interacting 
Bose gas. 
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FIG. 4: Structure factor hujS'/2fi (a) and TiwSi/2n (b) in the 
limit of small momenta k — > 0. The structure factor exhibits 
a peak at ui « 2/x, and is quenched above twice the bandwidth 
~ 24 J in the presence of an optical lattice (b). 

Subjecting the Bose gas to an optical lattice, we have 
to replace the free dispersion relation e k by the band 
structure of a particle in a periodic potential (in the tight 
binding limit, e k = 2J [3 — cos(fc K 6) — cos(k y b) — cos(k z b)} 
with the hopping energy J and the lattice constant b). 
In addition, Umklapp processes transferring momentum 
to the optical lattice have to be considered. Restricting 
the analysis to the lowest Bloch band and small momenta 
k — > 0, the expressions for the response function i|l(J|) and 
(|13|l remain valid and the contribution S^(lo, k) to the 



4 



structure factor takes the form 
2^ 2 



S£(w,k) = 



(huj) 4 J K (2ir 



^3 (k.Ve q ) 2 ^-2£ q ) (16) 



with K the first Brioulline zone. Furthermore, the 
structure factor turns quasi-periodic, S' L (u),'k + K) = 
\ok\ 2 S'^w , k) for k — > with K a reciprocal lattice 
vector and c K a numerical prefactor accounting for the 
shape of the wave function within a potential well (c.f., 
Ref. (l?! for a similar discussion of the quasi-periodicity) . 
In Fig. ^b) we show the structure factor S' l (Huj/2ii) = 
8m*n 2 S' lj (uj, k) / (n D ft 2 k 2 ) in the tight binding- and k — > 
limit with the effective mass m* = h 2 /(2Jb 2 ). 

In a recent experiment, Stoferle et al. |18| have deter- 
mined the energy transfer to an atomic gas via a modu- 
lation of the optical lattice with frequency u>. This mod- 
ulation perturbs the Bose gas with an external periodic 
potential of amplitude SV(u>,K) and K a reciprocal lat- 
tice vector. The energy per particle W transferred to the 
Bose gas during the time At takes the form (T = 0) 

W = ^ — 6V 2 (cj,K)nwS(cj,K). (17) 
2 n 

The main energy transfer observed in [lsj resides at large 
energies and is incompatible with the vanishing spectral 
weight in So at these frequencies in the following 
we attempt an explanation involving the quantum deple- 
tion. The structure factor S(u),K) has to account for 
the finite trapping potential: the wave functions describ- 
ing the Bogoliubov excitations involve momenta with a 
spread Ak. The structure factor in a finite trap derives 
from its counterpart in the homogeneous system via the 
average S(u>,K) w (S'(u,K + q))| q |<|Ak|- ln the inter- 
esting regime fku ~ fi, the structure factor S' involves 
two-particle excitations with large momenta and the fi- 
nite level spacing in the trap plays a minor role. In esti- 
mating the relevant spread Ak we have to account for the 
trap size R (Thomas Fermi radius) and the condensate 
healing length £ describing the rapid decay of the wave 
function at the trap boundary; making use of the exci- 
tations derived in Ref. [i(J , we obtain a spread of order 
|Ak| w l/(i?0 1/2 - Comb ining the quasi-periodicity of 
the structure factor in an optical lattice and the spread 
induced by the trapping potential, we find the estimate 



S(co,K)^\ Cli \ 2 (S[(iO, q )) H<lAkl 



(18) 



with S£(w,k) shown in Fig. Efb). In their experiments, 
Stoferle et al. [lsj report (among other) on the energy 
transfer to a trapped Bose gas subject to a 3D optical 
lattice residing in the superfluid phase. They find a pro- 
nounced peak in the excitation spectrum centered at a 
frequency Huj rs 2/i and with a width of order of the 
band width, in rough agreement with our result for S[. 
Accounting for the smearing in the structure factor due to 
the finite trap as discussed above, the transferred energy 



(|17fl assumes the correct order of magnitude. We then are 
tempted to associate the observations made in 0] with 
the quantum depletion of the Bose gas enhanced by the 
presence of the optical lattice. We note that the optical 
lattice applied in the experiments is already appreciable. 
Hence, we expect the weak coupling analysis above to 
reproduce correctly the qualitative features in the exper- 
iment, while a quantitative discussion requires to include 
higher order corrections in the parameter V na 3 . 

In conclusion, the structure factor of a weakly inter- 
acting Bose gas exhibits a characteristic behavior asso- 
ciated with the structure of the quantum depletion of 
the condensate. A measurement of the structure factor 
not only allows for the detection of the quantum deple- 
tion in atomic gases but provides additional non-trivial 
information on the ground state of strongly interacting 
Bose gases, thus opening up new possibilities in studying 
complex quantum phases in atomic gases. 

We thank T. Esslinger, M. Kohl, H. G. Katzgraber, 
S. Huber, and W. Zwerger for stimulating discussions. 
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